A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 

11 (11.03.1; 12.03.3; 12.04.1; 13.25.2 ) 



ASTRONOMY 

AND 
ASTROPHYSICS 



Regularity in the X-ray surface brightness profiles of galaxy 
clusters and the M—T relation. 



D.M. Neumann, M. Arnaud 

CEA/DSM/DAPNI A/Service d'Astrophysique, CEA-Saclay, 
L'Orme des Merisiers, Bat. 709, F-91191 Gif-sur-Yvette, France 
email: ddon@cea.fr , arnaud@hep.saclay.cea.fr 

Submitted January 5, 1999 



On 
On 

00 



> 
(N 
On 
O 
y—i 

O 

ON 
ON 



Oh' 
i 

o 



X 



Abstract. We used archival ROSAT observations to in- 
vestigate the X-ray surface brightness profiles of a sample 
of 26 clusters in the redshift range 0.04 < z < 0.06. For 
15 of these clusters accurate temperatures (kTx > 3.5 
keV) were available from the literature. The scaled emis- 
sion measure profiles look remarkably similar above ~ 0.2 
times the virial radius (?"vt20o)- O n the other hand a large 
scatter is observed in the cluster core properties. We fit- 
ted a /3-model (with and without excising the central 
part) to all the ROSAT profiles to quantify the struc- 
tural variations in the cluster population, unraveling a 
robust quadratic correlation between the core radius and 
the slope parameter (3. We quantified the shape of each gas 
density profile by the variation with radius of the logarith- 
mic slope, a n . The bi- weight dispersion of a n among the 
clusters is less than 20% for any given scaled radii above 
x = 0.2. There is a clear minimum spread at x = 0.3, 
which is related to the existence of a correlation between 
core radius and (3. These ensemble properties are insen- 
sitive to the exact treatment of a possible central excess 
when fitting the profiles. On the other hand the scatter is 
decreased when the radii are scaled to ryT200- 

The regularity we found in the gas profiles at x > 0.2 
supports the existence of an universal underlying dark 
matter profile, as already predicted by theoretical works. 
It suggests that non gravitational heating is negligible 
for clusters with temperature above ~ 3.5 keV. The very 
large scatter observed in the core properties favor scenario 
where Cooling Flows are periodically erased by merger 
events. Our results are consistent with the classical scaling 
relation between Mass and Temperature (M oc T 3 / 2 (1 + 
z)~ 3 / 2 ). Accordingly the spread in the reduced mass pro- 
files derived from the hydrostatic isothermal /3-model is 
small. 

Key words: Galaxies: clusters: general - Cosmology: ob- 
servations, dark matter - X-rays: general 
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1. Introduction 

There is growing observational evidence that the physical 
properties of galaxy clusters obey scaling relations. Evrard 
(1997) showed that most observed clusters of galaxies have 
a similar fraction of hot gas, as compared total mass, 
/gas ~ 0.060 /i~ 3 / 2 . This fraction sets a lower limit on 
the baryon fraction in these dark matter dominated ob- 
jects. Furthermore, studies focusing on the relationship 
between the X-ray luminosity of the hot gas (or intraclus- 
ter medium, hereafter ICM) and its X-ray temperature 
done by Markevitch (1998), Allen & Fabian (1998) and 
Arnaud & Evrard (1998) have revealed a relatively tight 
correlation between these two parameters. Recently Mohr 
& Evrard (1997) found that there exists a strong corre- 
lation between the X-ray isophotal radius and the X-ray 
temperature of the hot ICM. The study of Hjorth, Ouk- 
bir & van Kampen (1998), based on a sample of clusters 
with good X-ray temperature estimates and mass inferred 
from lensing observations, suggests that there is a tight 
mass-temperature relation. 

Such scaling laws are expected if clusters of galaxy 
form a homologous population and are an indication of an 
underlying structural regularity in the cluster population. 
Hence their study can give useful insight into the forma- 
tion and evolution of clusters and the underlying cosmo- 
logical scenario. In addition they can provide an efficient 
way to estimate some cluster properties, which are diffi- 
cult to measure directly. In particular the M-T relation is 
of special interest, in view of the cosmological importance 
of determining cluster masses and the possibility to mea- 
sure accurately the temperature with current and future 
X-ray telescopes. 

Structural regularity is expected on theoretical grounds 
(Teyssier, Chieze & Alimi 1997). Numerical simulations by 
Navarro, Frenk & White (1996, 1997) indicate that CDM- 
halos with masses spanning several orders of magnitudes 
follow a universal density profile, whose shape is indepen- 
dent of mass or cosmology. The corresponding density pro- 
file of the hot gas captured in these cluster halos can be 
reasonably well described by an isothermal /3-model , as 



2 



D.M. Neumann, M. Arnaud: Regularity in the cluster surface brightness profiles and the M-T relation 



shown by Eke, Navarro & Frenk (1998). The X-ray surface 
brightness profile reads, for such a model: 

,1 -3/3+1/2 



S(9) = So 



1 



which translates into the gas density distribution: 



n s (r)=n g0 1 + (r/r c ) 



-3/3/2 



(1) 



(2) 



The /3-model is known since the early Einstein observations 
(Jones & Forman 1984) to provide a good fit to the X-ray 
profile of galaxy clusters. However Makino, Sasaki & Suto 
(1998) recently found that the gas density profiles deduced 
from simulations are too steep to match with the typical 
core radii r c derived from these X-ray observations. More- 
over it is well known that there are large variations, from 
cluster to cluster, in the observed /3-model shape param- 
eters, (3 and r c . 

In this paper we want to check the regularity of the gas 
distribution in clusters, described by a /3-model . We se- 
lected for this study a sample of nearby clusters observed 
with ROSAT and we unravel a relation between the core 
radius and the slope parameter /3. We discuss the impli- 
cations of our results for the variations of the gas density 
profile shape and the existence of a mass-temperature re- 
lation. 

The paper is organized in the following way: we de- 
scribe in Sect. ||our sample of clusters of galaxies and the 
observations. Sect. |3] shows the surface brightness profiles 
of the clusters. In Sect. ^ we present our isothermal /3- 
fits. Sect. H shows the correlation between r c and f3. In 
Sect. ^ we describe the consequences of the correlation on 
the shape of the gas density profiles and quantify its vari- 
ations and in Sect. (?] we study the mass profile and the 
M-T relation. In Sect. || we discuss our results and give 
our conclusions. 

Throughout the paper we assume H = 50km/sec/Mpc, 
A = 0, Q = 1 (g = 0.5). 



2. X-ray observations 

We used X-ray imaging data retrieved from the ROSAT 
data base at MPE. The sample we built consists of clusters 
of galaxies found by Abell, Corwin & Olowin (1989) in the 
redshift range z=0.04-0.06, which were in the field of view 
of a ROSAT pointed observation with an off-axis angle 
less than 10 arcmin. The redshift range is small enough to 
avoid to look at major evolutionary effects within the sam- 
ple and the typical cluster size for these redshifts is well 
matched to the PSPC central field of view - mostly inside 
the rib structure of the instrument - or to the HRI field 
of view. Furthermore we only selected clusters whose ex- 
tended emission can be fitted with an isothermal /3-model, 
i.e. the cluster must show a clear center in the X-ray emis- 
sion. Finally the signal to noise ratio of the cluster must 
be high enough to allow accurate modeling of the data. 
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Fig. 1. The ROSAT X-ray surface brightness profiles of 
the clusters in our sample. The profiles are corrected for 
vignetting effects and background subtracted. Full line: 
clusters in the spectroscopic subsample. 



The list of the 26 clusters selected is given in Tab.[j], 
as well as the exposure times and the detector used. The 
sample covers a large variety of morphological types, from 
relaxed spherical symmetric clusters like A2107, to clus- 
ters with substructures like A754 or A3559 and contains 
non-cooling flow clusters (e.g. A119) as as well clusters 
known to have strong cooling flows (e.g. A85). 

For the PSPC data we only used photons in the band 
0.5-2.0 keV, for the HRI we only took into account channel 
2-9, in order to optimise the signal-to-noise ratio. For each 
data set we calculated the exposure map, using the soft- 
ware implemented in the X-ray software package EXSAS 
developed at MPE for the PSPC data and the software 
developed by S. Snowden (Snowden 1998) for the HRI 
data. 

We searched the literature for the emission-weighted 
mean temperature of the clusters in our sample. We only 
considered temperature obtained with wide energy band 
experiments (Einstein/MPC, EXOSAT, GINGA or ASCA) 
to avoid systematic errors and exclude temperature esti- 
mates with statistical errors greater than 50%. The adopted 
temperature values and corresponding references are given 
in Tab. [l] for 15 clusters, which constitute our spectro- 
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Fig. 2. The derived emission measure profile of the clus- 
ters in the spectroscopic sub sample. The dotted line 
shows, for comparison, the emission measure profile of 
Abell 2163 (z=0.201, kT=14.6 keV, Elbaz, Arnaud & 
Bohringer 1995) 
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Fig. 3. The scaled emission measure profile (Eq. ^) of 
the clusters in the spectroscopic sub sample. The radius is 
normalised to ryT200 (Eq. ||). Dotted line: Abell 2163 pro- 
file. Beyond r/ryT200 ~ 0.1 the profiles look remarkably 
similar. 



scopic sub-sample. Note that this sub-sample does not 
contain low temperature clusters (kT > 3.5keV), because 
only the brightest clusters of our sample have good tem- 
perature measurements (see Fig. |l|). All these clusters have 
flux greater than 1.7 10~ n ergs/s/cm 2 . The spectroscopic 
sample is 80% complete at that flux limit, when compared 
to the complete sample derived by Ebeling et al. (1996) 
from the ROSAT All-Sky Survey. 



3. X-ray surface brightness profiles 

Fig.^ shows the vignetting corrected X-ray surface bright- 
ness profiles of all the clusters in the sample. We bin the 
photons into concentric annuli centered on the maximum 
of the X-ray emission. For the PSPC we use a width of 15 
arcseconds per annulus and a total of 200 annuli. For the 
HRI we use a width of 10 arcseconds per annulus and a 
total of only 100 annuli, due to the smaller field of view of 
the HRI. We cut out serendipitous sources in the field of 
view or cluster substructures, if they show up as a local 
maximum. The background was subtracted using data in 
the outer part of the field of view. 



The emission measure along the line of sight at radius 
r, EM(r), can be deduced from the X-ray surface bright- 
ness, S(6): 

4tt (1 + z) 4 S(9) 



EM(r) 



d A (z) 



(3) 



A(T,z) 

where A(T, z) is the emissivity in the ROSAT band, tak- 
ing into account the interstellar absorption and the instru- 
ment spectral response, and d\{z) is the angular distance 
at redshift z. A(T, z) depends, although weakly, on the 
cluster temperature and redshift. The emission measure 
is linked to the gas density n g by: 

n 2 JR) RdR 



EM(r) 



(4) 



The shape of the surface brightness profile is thus gov- 
erned by the form of the gas distribution, whereas its nor- 
malization depends also on the cluster overall gas con- 
tent. If clusters formed a structurally similar population, 
all the brightness profiles would appear as parallel curves 
on Figure [l] (they would differ by a translation in the log- 
log plane). More precisely similarity means that the clus- 
ters constitute a one parameter population: each cluster 
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Table 1. ROSAT observation analysis summary, p/h stands for PSPC or HRI. All quoted errors are 2 — a errors. The 
redshifts are taken from NED. When the best fit /3-model A and B are identical, we only give the parameters once. 
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: Markevitch 1998; a : Arnaud & Evrard 1998; d : David et al. 1993; ^ : Yamashita (priv. comm.); : Edge & Stewart (priv. 



comm. quoted in White, Jones & Forman 1997) 



is characterized by its physical dimension a^, and the dis- 
tribution of any given physical quantity Q is described 
by a dimensionless function Q, common to all clusters: 
Qi{f) = Q(o-i)Q{x), where x is the scaled radius x = rjai. 
The normalization factors Q(ai) define the scaling rela- 
tions between the global quantities (e.g mass, mean tem- 
perature, luminosity . . . ) in the cluster population. 

We first examined if our data were consistent with the 
similarity expected in the simplest model. In the spherical 
collapse model, the overall density contrast S of virialized 
objects is fixed and is of the order of 200 for £7 = 1, and 
the natural scaling radius is the virial radius defined as the 
radius containing this density contrast^]. Here the density 
contrast is S = p/p c , where p c {z) is the current critical 



density of the Universe, p c (z) = 3H /8"7rG(l- 



and p 



is the mean cluster density. The virial radius is not directly 
measurable. However, if, in addition, one assumes struc- 
tural similarity, the virial theorem provides a scaling rela- 



1 Note that numerical experiments support that definition of 
the virial radius: the edge of the virialized part of clusters is 
indeed found at a density contrast ~ 200 



tion between the virial mass (Mvt20o) and radius (rvT20o) 
and the overall X-ray temperature Tx: AfvT200 / r VT200 °c 
Tk, whereas by definition MvT20o/(4/37r,o c (z)ry T200 ) = 
200. This leads to the well known scaling relations: 

r V T200 « (l + zr 3/2 T x /2 (5) 

Af V T200 oc (1 + z)- 3 / 2 T 3/2 (6) 
The emission measure (Eq. ||) then scales as 

EM{r) oc / 2 T x /2 (1 + zf/*EAl{x) (7) 



where x — r/rvT200 is the scaled radius, and EM is a 
dimensionless function, the same for all clusters. The gas 
mass fraction / gas is not necessarily a constant but sim- 
ilarity implies that it only depends on cluster mass (or 
equivalently temperature). 

We thus considered the observed scaled emission mea- 
sure profiles EMx(x): 

S(0(x)) 



EM x (x) = 



(l + z)V2 Tx /2 A(T) 



0{x) 



fVT200 / x 

x , , , (8) 



d A (z) 



defined from Eq. and Eq. ^ ignoring possible variation 
of /gas with Tx- To compute the scaling radius ryT200 we 
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Table 2. Description of the /3- model 



/3-model description 



A isothermal /3-model fitted to the entire 

cluster surface brightness profile 
B isothermal /3-model excising the central region 

when an improvement of \ 2 can be achieved 



used the normalisation factor obtained by Evrard, Metzler 
& Navarro (1996) from numerical simulations: 

1/2 

rVT2oo = 3690 ( T7 f^- ) (1 + z)^' 2 kpc (9) 



T x 
10 kev. 

Note that this normalization is simply for convenience, the 
exact calibration of the scaling relations does not matter 
to check similarity of the profiles. 

The scaled emission measures profiles EMx(x) are shown 
on Fig. U for the clusters in the spectroscopic sub-sample. 
They can be compared to the corresponding unsealed pro- 
files EMx(r) plotted on Fig. |^. The scaled profiles are 
clearly not identical, but several features are striking. First, 
we detect the emission nearly up to the virial radius. We 
are effectively studying the large scale density distribu- 
tion of the gas, and not just the central cluster core. Sec- 
ondjjwhile the dispersion is very large at low scaled radii 
(EMx at the center span nearly 2 orders of magnitude), it 
rapidly decreases with radius. Beyond r/rvT200 ~ 0.1 the 
profiles look remarkably similar and agree to within about 
a factor of two. We also note that the scaling procedure 
has significantly decreased the difference between the pro- 
files, but only in that external region (compare Fig. || and 
Fig. H). The unsealed emission measures show a relative 
standard deviation of about 96% at 950 kpc and 75% at 
1500 kpc. In comparison, the relative standard deviation 
of the scaled profiles is 41% at r/ryT200 = 0.3 and 39% 
at r/rvT200 = 0.5, i.e the scaling has reduced the scatter 
by about a factor of two. As the temperature range of our 
spectroscopic sub-sample is limited (3.5 < kT < 9.1keV), 
we plotted for comparison the scaled profile of A2163. It 
is one of the hottest cluster known (Arnaud et al. 1992) 
and is located at higher redshift than the clusters of our 
sample. It fits remarkably well within our sample. 

In summary, the scaled emission measure profiles show 
indications for similarities in outer regions (r > 0.2 tvt2qo)- 
On the other hand in the centers the profiles show a large 
dispersion, which we will discuss later on. To further quan- 
tify the density profile shape and investigate its possible 
variations within the cluster population, we fit, in what 
follows, an isothermal /3-model to each cluster profile. 



4. Fitting the isothermal /3-model 

We performed two different /3-modcl fits for each cluster. 
First we determined the /3-model that fits at best the en- 
tire cluster emission (hereafter /3-model A). However it is 
well known that the overall /3-model is a poor description 



of the central region of some clusters where excess emis- 
sion is observed (due to a cooling flow, a central point 
source . . . ) . We indeed got in several cases very large x 2 
values (see Tab [|). Hence, we also tried to minimise the 
reduced x 2 by excluding the central bins from the fit. The 
best fit /3-model (hereafter /3-model B) was determined 
by increasing the radius of the cut-out region, i? cu t , up to 
the radius where the reduced x 2 stopped decreasing. We 
limited i? C ut to 350 kpc, as a too large excluded region in 
the center precludes the determination of the core radius. 
Moreover it is safe to assume that the contribution of any 
Cooling Flow excess is negligible beyond such radius. 

A brief description of /3-model A and B is given in 
Tab 0. The best fit parameters, (3 and r c , are listed in 
Tab [j] for both models, together with the i? cu t values. For 
about 1/4 of the clusters, excluding central bins does not 
improve the fit. In that case the best fit /3-models A and 
B are identical and we only give the parameters once in 
the Table. 

For the other clusters, excluding the central part yields 
on average larger /3 and core radii values (see Tab. Q). This 
is a well known effect. If one tries to fit a /3-model to the 
entire profile when there is a central excess, too small core 
radius and /3 values are derived: decreasing r c allows us to 
fit better the central part, while decreasing /3 compensates 
for the subsequent flux deficit induced at large radii. Al- 
though the /3-model B is better in that respect, it has also 
its drawback. If the excluded region is larger than strictly 
necessary to avoid the central excess, the determination of 
the core radius is degraded: the uncertainty is increased 
and the best fit value can be biased towards large values. 
We will thus consider both models in the following. They 
can be viewed as two extreme /3-models , allowing to as- 
sess the impact of the exact treatment of the cluster core. 
This is important since the core properties vary greatly 
from cluster to cluster, as shown in the previous section. 

For A754, which is undergoing a merger phase, as indi- 
cated by many authors (Zabludoff & Zaritsky 1995; Hen- 
riksen & Markevitch 1996; Roettiger, Stone & Mushotzky 
1998), we varied the center used for the concentric binning 
and adopted the profile that provides for the global fit the 
lowest reduced x 2 . 

The most compact clusters A780 and A1991 were ob- 
served with the PSPC, and their core image is somewhat 
blurred by the PSF (« 20-30 arcsec FWHM). For A780 we 
estimated that the core radius is overestimated by about 
10%, and for A1991 by around 20%. For the PSPC point- 
ing of A85 (the cluster which has after A780 and A1991 the 
third smallest angular core radius) the effect of the PSF 
is of the same order as the statistical uncertainties. For all 
other clusters, which either have a larger core radius or 
were observed with the ROSAT HRI Q we calculated that 



2 the HRI has a PSF with a FWHM of 4-5 arcseconds, which 
corresponds to about 4-8 kpc 
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Fig. 4. Slope parameter [3 versus core radius r c for the 26 
clusters of the sample (Tab.|l|). The isothermal /3-model 
has been fitted to the entire ROSAT surface brightness 
profile of each cluster (/3-model A). The errors bars are 
2-cr errors. /3 and r c are correlated. The full line shows 



1 + (r c /r B Y 



with 



the best fit parabolic relation: [3 — (3ci 

Po = 0.55 and r s — 885 kpc. The dotted line shows the 
parabolic correlation obtained for the /3-model B param- 
eters: Po — 0.61 and r s = 1020 kpc (see also Fig. ||). 

the effects of the PSF are much smaller than the statistical 
errors. 

All clusters with large error bars, such as A76, A3093, 
A3223 and A3559, have clearly visible substructures. 



5. Correlation between the core radius and (3 

5.1. The (3 and r c values and the significance of the 
correlation 




500 



1000 



r in [kpc] 



Fig. 5. Same as Fig. [| only that we now display the fit 
results of /3-model B: isothermal /3-model excising the 
central parts of the profile, if an improvement of reduced 
X 2 can be achieved (see Tab.|l|). Clusters belonging to the 
spectroscopic subsample are marked with filled circles. 



alent plot for /3-model B. A striking feature is that the 
data points are not distributed all over the r c -/3 plane. 
There is a clear trend of increasing slope parameter with 
core radius, although the correlation is not very tight, spe- 
cially at low core radius values. A quantitative measure of 
the correlation significance was obtained from three dif- 
ferent tests, the Pearson's test, the Spearman's test, and 
the Kendall's r test (see Press et al. 1993). All three tests 
(applied to fit A results) give a high correlation signifi- 
cance, higher than 99.97%. The highest correlation signif- 
icance is obtained with the Pearson's test, with a value 
of 1. — 2.99 x 10 -10 . The corresponding Pearson's rank 
correlation is r = 0.902. 



The slope parameter /3 is not the same for all clusters, as 
would be the case for perfectly similar /3-model density 
profiles. The (3 values span a wide range, 0.47 — 1.04 and 
0.52 — 1.4 for the /3-model A and B, respectively. The 
distribution is not symmetric, it is clearly concentrated at 
low values, with rarer high values: the median value is 0.59 
for A (0.68 for B), ~ 85% of the values fall in the narrower 
range 0.47(0.52) — 0.75 and the standard deviation is only 
0.14(0.21) or 24%(30%) of the median value for model A 
(model B). The same general features are observed for the 
core radius distribution. The median values are 143 kpc 
and 278 kpc for models A and B, respectively. 

Fig. ^ displays the slope parameters (3 against core ra- 
dius r c , obtained with the /3-model A. Fig. wis the equiv- 



5.2. The form of the r c -(3 relation 

The highest score obtained with the Pearson test indi- 
cates that the correlation between r c and /3 is close to 
being linear. One notes that the slope parameter is nearly 
constant at low core radii and starts to increase with r c 
above r c ~ 250kpc only, which suggests a parabolic rela- 
tion between the two quantities. 

We examined which relation, linear or parabolic, bet- 
ter accounts for the data. For conciseness, the following 
discussion is based on a quantitative analysis of the pa- 
rameters of the /3-model A (similar results are obtained 
with /3-model B). 



1.4 
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Table 3. The best fit parameters of the /3-core radius 
relations 




0.2 0.3 

r c / r VT200 

Fig. 6. The correlation between scaled core radius x c = 
r cAvT200 and (3 for the 15 clusters of the spectroscopic 
subsample (see also Tab.0). The data points are for the 
/3-model B and the error bars take into account the un- 
certainties in rvT200- Full line: (3-x c relation for /3-model 
B. Dotted line: (3-x c relation for /3-model A. We also dis- 
play the results of A2163 located at z = 0.201 (Arnaud et 
al. 1992; Elbaz, Arnaud & Bohringer 1995) and CL0016 
at z = 0.54 (Hughes & Birkinshaw 1998). 



(10) 



We considered the two parametric functions 
= P l o [l + (r c /rj)] 
and 

/3 = /3 [l + (r c /r s ) 2 ] (11) 

and fit them to the data. Since both quantities (3 and r c 
are marred with errors and the uncertainties are highly 
correlated, there is no well established method to define 
the best fit function parameters and to quantify the good- 
ness of the fit. We used the following empirical least- 
squares method. Best fit parameters were derived by omit- 
ting the errors on j3 and r c . The x 2 was minimised using 
the downhill simplex method (see Press et al. 1993). We 
got (3q = 0.49 and r l s — 870 kpc for the linear function 
and (3q = 0.55 and r s = 885 kpc for the parabolic function 
(the corresponding curve is plotted on Fig.||) . To compare 
the goodness of the fits obtained with the two functions 
we performed a subsequent x 2_ test including the error on 
the j3 quantity only. We utilised the 2 — a errors instead of 
the 1 — a errors. As the uncertainties of (3 and r c are cor- 
related we assumed this could crudely take into account 
the two errors. 

The reduced % 2 values for the entire cluster sample are 
14.6 and 9.2 for the linear and parabolic function respec- 
tively. The fit is formally better for the parabolic relation 



Relation 


/3-model 


normalisation 


scale 


/8-Tc 


A 


0o = 0.55 


r s = 885 kpc 


/9-Tc 


B 


/3o = 0.61 


r s = 1020 kpc 


(3-Xc 


A 


I3' = 0.53 


:r s = 0.27 


f3-x c 


B 


Po = 0.59 


x a = 0.34 



but both x 2 are very high. If we exclude four clusters, A85, 
A780, A2319, and A3562, the reduced X 2 drops to 1.68 for 
the parabolic function. These four clusters lie offset of the 
correlation and have extremely small statistical error on 
(3, less or equal to 2%. They boost therefore the x 2 value. 
The /3-model parameters of these four clusters change sig- 
nificantly when one cuts out the central part (see TabJIJ) 
and it is likely that using the global fit for these clusters 
is not a good approach. On the opposite we still obtain a 
reduced x 2 of 3.36 with the linear relation after the ex- 
clusion of the four principal outliers of this relation. Only 
when excluding seven clusters, we do obtain a reduced x 2 
of 1.78, which is still higher than the x 2 of the parabolic 
fit, with four clusters excluded. 

While keeping in mind that our criteria are not rigor- 
ous, we conclude that a parabolic relation between r c and 
(3 is a better description of the correlation between these 
two quantities than a linear relation. Henceforth, we will 
only consider a parabolic relation from now on. 



5.3. The best fit r c -(3 relation 

We determined the best fit parabolic relation r c -(3 for both 
/3-models A and B. The robustness of the r c -(3 relation 
must be noted; it is not sensitive to the exact treatment of 
the central part of the cluster profile. From the /3-model 
B data, plotted on Fig. g, we got /3q = 0.61 and r s = 
1020 kpc. These values are close to those obtained for the 
/3-model A: (3 is increased by only ~ 10% and r s by ~ 
14%. The reduced x 2 is smaller: x 2 = 1-92 for the complete 
cluster sample; the x 2 drops to 1.42 if one excludes A780. 
Part of the improvement is certainly an artifact due to 
the larger uncertainties on (3. The two relations look very 
much alike, as can be seen on Fig. |^. 

As (3 is a dimensionless quantity, while r c is not, one 
would rather expect a relation between (3 and a scaled 
core radius. Fig. || shows the correlation between [3 and 
the scaled core radius x c — rc/rvT200> where the virial 
radius ryT200 is given by Eqj9| and the data points are 
from the /3-model B. We fitted the parabolic relation 



f3 = (3' l + (x c /x s f 



(12) 



to the cluster parameters of the 15 clusters in the spectro- 
scopic subsample. We obtained /3q = 0.59 and x s = 0.34 



for the /3-model B and similar results (3' 
x s = 0.27 for the /3-model A (see also Tab|| 



0.53 and 



8 



D.M. Neumann, M. Arnaud: Regularity in the cluster surface brightness profiles and the M-T relation 



Table 4. The (3-r c relation for various cluster samples 

Number of Choice /3-model /3o r a 

clusters criteria 

~~26 none A 055 885 kpc 

26 none B 0.61 1020 kpc 

~ 13 r det > 1.45Mpc A 6\54 859 kpc 

13 x 2 < 1-6 A 0.56 980 kpc 

13 r det > 1.45Mpc B 0.60 1050 kpc 

15 x 2 < 1-5 13 0.64 1080 kpc 



the quality of the data and the cluster region used for the 
fit. We have already mentioned that a central excess can 
yield artificially low (3 and core radii. At large radii the 
/3-model tends to a simple power law. If instead the den- 
sity profile continuously steepens, the /3 value obtained 
from the fit will increase with the outermost radius used 
in the fit (see for instance the numerical simulations of 
Navarro, Frenk & White 1995) and as a result probably 
also the derived core radius. Moreover the magnitude of 
these effects will depend on the accuracy of the /3-model 
approximation, which is likely to vary from cluster to clus- 
ter. All together this could induce a spurious correlation 
between the derived parameters, similar to the one ob- 
served, i.e. an apparent increase of (3 with r c . 

To see if this a serious issue we first investigated how 
the f3-r c relation is affected by the size of the region used 
for the fit. We can first compare /3-models A and B. We 
also derived the relation for the half sample of the 13 clus- 
ters with the largest rdet, where rdet is the maximum ra- 
dius at which the X-ray emission is detected with a sig- 
nificance greater than 3o\ Second we investigated if the 
relation is sensitive to the quality of the fit obtained with 
the /3-model . For this purpose we selected the clusters 
with the lowest reduced % 2 and derived the corresponding 
0-r c relation. The \ 2 value is an indicator of the accuracy 
of the /3-model approximation (but also depends on the 
quality of the data). 

The parameters of the relations derived for the vari- 
ous cluster sub-samples are given in Tab. |j. In all cases, 
a parabolic (3-r c relation well accounts for the data. Fur- 
thermore for each model (A or B), the parameters /3o and 
r s depend very weakly on the selection on rdet (0.3 — 3% 
effect) or on the selection on the quality of the fit (2 — 10% 
effect). The largest discrepancy is actually seen between 
/3-models A and B, i.e if one includes or not the central 
excess region in the fit. The effect is small as discussed in 
the previous section. From this observed robustness of the 
(3-r c relation we concluded that it is unlikely that it is an 
artifact due to the imperfection of the /3-model . 



5.4.2. The fitting process 

To check the validity of the fitting procedure, we con- 
struct artificial clusters, which lie outside the correlation 
and have exposure times and countrates similar to the 
ones in our sample. We added Poisson noise to the im- 
ages and treated these artificial clusters in the same way 
as our real data and fit the isothermal /3-model to them. 
The intrinsic values for [3 and r c always lie well within 
the 2 — (T-errors of the fit results of the simulated clusters 
images. Thus, the fit procedure is correct and could not 
introduce a spurious correlation. 



An important question is whether (3 is more tightly 
related to r c or to x c . We cannot definitely answer this 
question with the present spectroscopic subsample. It is 
too small and covers a too narrow range in temperature 
and redshift. There is only a factor 2.4 between the mini- 
mum and maximum temperature. As rvT200 only scales as 
-\/Tx, the effect of the scaling cannot be dramatic. How- 
ever there are several indications that the 'scaled' relation 
is the most relevant one, as expected. First we compared 
Fig. ^ to Fig. H at large (3 values (/3 > 0.65), where the de- 
pendence of (3 with r c is maximal. For the clusters in the 
spectroscopic subsample (filled circles) the scatter around 
the best fit relation has been decreased by the scaling pro- 
cess. Second we considered test clusters outside the T and 
z range of our sample. A2163, the highest temperature 
cluster, fits well with both the (3-x c relation and the /3- 
r c relation, so this is not conclusive. On the other hand 
C10016+16, the most distant cluster with good estimates 
of (3 and r c ((3 = 0.728, r c = 298 kpc) fits much better 
with the f3-x c (see Fig. ||). This issue is further discussed 
in Sect. ||. 

A summary of all fits of the (3-r c relation is given in 
Tab. |3|. We do not give errors on the relation parameters, 
/3o and x s (or r s ). We do not know of any proper method 
to estimate them since the uncertainties on the related 
quantities (3 and r c are highly correlated. We only suggest 
that they could be of the order of the difference observed 
when considering the two /3-model A and B. We also want 
to stress that the relation parameters we determine here 
arc in very good agreement with the direct study of the 
density profile shape presented below in Sect. ^. 

5.4- Checking on spurious effects 

Could the derived correlation be an artifact of our data 
analysis? This question must be addressed since our re- 
sults could in principle be affected by i) a choice of an 
inadequate model ii) the fitting process iii) systematic er- 
rors in the modeling of the instrument characteristics. We 
now discuss each point in turn. 

5.4.1. The model 

The /3-model is only an approximation of the real gas den- 
sity profile. Hence, the parameters derived may depend on 
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5.4.3. The instrument 

As already discussed above, only two clusters could be 
affected by the finite instrument PSF, namely A780 and 
A1991. For the others, the PSF effects are of the same or- 
der than the statistical uncertainties or smaller. Excluding 
these two clusters from the fit changes the results of /3q and 
r s by less than 3%. 

We also verified the exposure map calculation (or vi- 
gnetting effects). We used the observations of clusters, 
that were not included in our sample because we found no 
prominent extended cluster emission. This is the case for 
A195 (observing title MKN 359) observed by the PSPC, 
and A1213 observed by the HRI (observing title 1113+29). 
For both pointings we divided the image of the actual data 
by the exposure map. There is no significant variation in 
the overall photon distribution when excluding obvious 
sources: the countrate is flat over the whole field of view 
of the corresponding detector. This also validates our im- 
plicit assumption of a flat background in the FOV. 



6. The shape of the gas density profiles 

A perfect similarity between the density profiles is ruled 
out by the variation of [3 from cluster to cluster. On the 
other hand the correlation found between the two shape 
parameters of the density profile indicates some structural 
regularity in the gas distribution, that we study in this sec- 
tion. For that purpose we use the density profiles derived 
from our /3-model fits of the observed brightness surface 
profiles. 



6.1. Quantifying the shape: the slope profile 

The shape of the gas density distribution may be quanti- 
fied by the variation with radius of its logarithmic deriva- 
tive. This derivative is also important for the estimate of 
the mass profiles (see Sect. 0). For the /3-model form, it 
reads: 



a n (r) 



d\ogn g (r) 
dlogr 



3 (3 



1 + (r c /rY 



(13) 



or equivalently, in term of the scaled radial coordinates 

3 (3 



a n (x) 



1 + (x c /xY 



(14) 



We will refer to this function as the slope profile: a n (r) 
is the slope of the gas distribution in the log-log plane at 
radius r. It depends on the two shape parameters of the 
/3-model : the core radius and /3. 



6.2. The slope profiles for a theoretical parabolic (3 - core 
radius relation 

For a parabolic relation between (3 and r c (Eq. |ll[), the 
slope profile may be written as: 



a n (r) = 3 (3 



1 + (r c /r B y 



1 + (r c /rY 



(15) 



Similarly for a parabolic relation between (3 and the scaled 
core radius: 



a n (x) = 3 (3' 



i + (x c /x s y 



1 + (x c /x)" 



(16) 



The physical implications of the parabolic relation, is 
straightforward from these two equations: i) if (3 is related 
to r c , there is a specific radius, r s , where the slope of the 
gas density profile is the same for all clusters, ii) if (3 is 
instead related to the scaled core radius, this occurs at a 
specific scaled radius, x s . 

Let us first assume a relation between [3 and r c . At 
the specific radius r s the slope is equal to 3 [3q, where 
[3q is the normalisation factor of the parabolic relation. 
At very large radii, r >> r c , a n (r) tends to 3 j3 - the 
/3-model is equivalent to a power law - and there is no 
impact of the f3—r c relation. At low radii, a n (r) roughly 
scales as (r c /r) and in view of the large dynamical range 
of core radii values, it can vary by more than an order of 
magnitude from cluster to cluster. This is illustrated on 
Fig. ^ where we have plotted with thick lines two extreme 
slope profiles corresponding to the best fit f3— r c parabolic 
relation derived in Sect. (Eq. |l5| with parameters from 

Tab. ||). The first one, (r), corresponds to the min- 
imum allowed [3 value: (3 = /3q or r c = 0. In that case 
a n (r) keeps constant with radius: a[^\r) = 3 [3q. The 

(2) 

second one, a n (?")> corresponds to the maximum (3 value 
observed in our sample. The two curves cross at r s . 

The same properties hold for the (3-x c relation, in 
terms of the scaled radial coordinates. 



6.3. The observed slope profiles 

The slope profiles a n (r) derived from model A and B are 
plotted on the left panels of Fig. ^ for all clusters in our 
sample. On the right panels are plotted the profiles a n (x) 
against scaled radius for the spectroscopic subsample. 

To quantify structural variations within the cluster 
sample, we estimated at each radius, the center and spread 
of the distribution of slopes at that radius. We used both 
the classical estimators (mean value and standard devia- 
tion) and the bi-weight estimators (bi-weight location and 
scale). The latter are more resistant and robust than the 
former (Beers, Flynn & Gebhardt 1990). The results are 
shown in Fig. M. The "average" slope (top panels) and 
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Fig. 7. The variation with radius of the logarithmic slope of the gas density distribution. Top-left panel: the slope 
profiles a n (r) = —d\og(n g (r)/d\ogr versus radius r derived from the best-fit /3-model A, for all clusters of the sample. 
The bold lines correspond to the extreme slope variations allowed by the j3— r c correlation (Eqjll] with parameters from 
Tab. H and Eq.[l5]): the horizontal line is for /3 = fio = 0.55, r c — and the other line is for /3 = 1.4, r c = 1100 kpc. 
Top-right panel: the slope profiles a n (x) = —d\og(n g (x)/d\ogx versus scaled radius x = r/rvT200 for the spectroscopic 
subsample (/3-model A). Bottom panels: Same as top panels but for /3-model B. In that case the profiles are plotted 
only beyond the excised central region. 



the "relative dispersion around the mean" , defined as the 
spread normalised to the average value (bottom panels) 
are plotted as a function of radius (left) for the whole 
sample and scaled radius (right) for the spectroscopic sub- 
sample. The results of model A (dotted lines) and model B 
(full lines) are provided. We also computed the bi-weight 
means of the core radii and cut-out radii R cu t, the loca- 
tions of which are indicated by arrows in the figure. 

6.4- General properties 

The characteristics of the observed profiles are those ex- 
pected from the best fit /3-core radius relations we derived. 
Let us consider for instance the slope profiles a n (r) derived 
from model B. At small radii the observed profiles lie well 



within the two extreme curves ai (r) and a4 (r) corre- 
sponding to the best fit j3-r c relation and the scatter is 
very large among the cluster population (Fig. |?], bottom- 
left panel). This scatter diminishes with radius until it 
reaches a minimum at the specific radius r s — 1020 kpc 
(Fig. [7| and Fig. || bottom-left panel). At that radius the 
average a n value agrees well with the best fit 3 (3q value 
(Fig. H top-left panel). The classical and bi-weight estima- 
tors give the same results, the slope distribution at that 
radius being close to Gaussian. The residual scatter on a n 
at r s - the slope is not exactly the same for all clusters - 
can naturally be explained by the scatter around the best 
fit j3—r c relation. Above r s the scatter increases again, but 
saturate at a value of the same order as the variations of 
3 P, the asymptotic value of the slope. For instance, at 2.5 
Mpc the bi-weight mean slope is 1.96 and the bi-weight 
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Fig. 8. Structural variations within the cluster sample. The average of a n (top panels) and its relative dispersion 
(bottom panels) are plotted as a function of radius (left) for the whole sample and scaled radius (right) for the 
spectroscopic subsample. The relative dispersion is defined as spread normalised to the average value. Thick lines: bi- 
weight estimators. Thin lines: classical estimators. Dotted lines: results for /3-model A. Full lines: results for /3-model 
B. Arrows in left-top panel: They indicate the bi- weight mean values of the core radii r Ci A and r Cj B for /3-model A and 
B respectively; the bi-weight mean i? cu t; the specific radius r Si A(r S) B) of the f3-r c correlation for /3-model A(B). Arrows 
in right-top panel: Same for scaled coordinates. Dashed line in top-right panel: slope profile obtained for (3 = 0.67, 
x c = 0.13. 



relative spread is 17%, very close to the corresponding val- 
ues for 3 (3: 1.98 and 18% respectively. At these large radii 
the bi-weight mean and dispersion are smaller than the 
classical estimators, reflecting the non-Gaussian (3 distri- 



bution which presents a tail at high values (see Sect 5.1) 



These general properties hold both for /3-model A and 
B. They do not depend on the choice of the radial coordi- 
nates (physical or scaled radius). However, there are some 
quantitative differences between the corresponding data 
sets. 

The relative dispersion of the slope a n (r) is higher for 
the /3-model B than for the /3-model A at all radii. This 
may result in part from the larger uncertainties on (3 and 
r c , which introduce extra scatter. Using the scaled radii 



instead (Fig. || bottom-right panel), the relative disper- 
sion of a n for /3-model B decreases. The minimum dis- 
persion, about 9%, is similar for A and B in this case. 
This further supports the idea that it is better to work 
with scaled coordinates, in agreement with the results of 



Sect. 5.3. Furthermore the bi-weight dispersion for model 



B now becomes slightly smaller than the dispersion for the 
model A, beyond the region where the dispersion is mini- 
mal (x ~ 0.3). This is more satisfactory since the model B 
should be a better fit to the data at large radii as model 
A is expected to introduce extra scatter. From now on we 
will only work in scaled coordinates. 

The mean slope value is smaller at low radii (the den- 
sity distribution is flatter on average) and higher at high 
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radii (steeper density distribution) for the /3-model B than 
for the /3-model A (Fig. || top panels) . This is a direct con- 
sequence of the treatment of a possible central excess, as 
explained in details in Sect. [|. This central excess typi- 
cally occurs within x cu t = 0.08, the bi-weight mean cut- 
out radius. Beyond r/rvT200 ~ 0.1 the systematic differ- 
ence between the two models becomes smaller than the in- 
trinsic dispersion observed within the cluster population. 
The mean slopes for model A and B became identical in 
the region of minimum dispersion (x ~ 0.3). This is con- 
sistent with the robustness of the /3— x c relation already 
noted (Sect. 5. 3). This means that the determination of 



the mean slope around that scaled radius is very robust 
and does not depend on the exact model used to fit the 
data. 

The mean slope profile can be approximated within 
a few percent with the profile corresponding to the bi- 
weight mean core radius x c = 0.13 and /3 value j3 — 0.67 
(see Fig. |§ top- right panel). 



6. 5. The standard density profile and structural 
variations in the cluster population 




r/r 

VT200 

Fig. 9. The normalised mass profiles of the 15 clusters in 
the spectroscopic subsample. For each cluster the mass 
is derived from the hydrostatic equation and our best fit 
isothermal /3-model (Dotted line: /3-model A; Full line: 
/3-model B). The mass is normalised to Mvt2Q0) the mass 
at density contrast 200 derived from the numerical simu- 
lations of EMN96. The bold line is the profile derived by 
these authors. 



In view of the shape characteristics of the observed density 
profiles, it is natural to define a standard profile: this den- 
sity profile follows a /3-model in scaled coordinates with 
x c = 0.13 and /3 = 0.67. We recall that the cluster physical 



radius scales as T- 



1/2 



The typical deviations from this standard profile are 
the following: 

— A central excess compared to the /3-model can be ob- 
served in the central part, typically within x cut = 0.08. 
In that region the scatter in density profile shape is 
very large (more than 100% variation in logarithmic 
slope) among the cluster population. 

— Beyond x = 0.2 the bi-weight dispersion of the ob- 
served slopes a n (x) becomes less than 20% of the slope 
of the standard /3-model profile. 

— The dispersion is minimal at x = 0.3 where a n (x) — 
1.67 ± 0.18. This corresponds to the observed correla- 
tion between core radius and (3 values. 

Note that these values are derived from our spectroscopic 
subsample (Fig. ||), which concerns relatively hot clusters. 
Cooler clusters may present larger variations. 



7. Consequences for the total mass profile and 
the M-T relation 

X-ray astronomers usually utilize the hydrostatic approach 
to calculate the total mass of a cluster within a certain ra- 
dius 0: 



If the ICM is isothermal or if the temperature gradient 
can be neglected in front of the density gradient, one may 
express the total mass as a function of the X-ray temper- 
ature and the gas density logarithmic slope a n . Using the 
scaled radial coordinates one gets: 



M tot (x) = 4.13 x 10 



Mtot(r) = - 



YT fd\ogn g {r) dlogT 
G/inip \ dlogr dlogr 



(17) 



15 



3/2 / n 

x a n (x) 



M Q (18) 



lOkeVy (l + z)3/ 2 

where a n (x) is given by Eq. |l4] for a /3-model . 

It is convenient to normalize the mass profile to the 
virial mass at density contrast 200 derived from numerical 
simulation (Evrard, Metzler & Navarro 1996): 

3/2 

M\ , -mm : 2.92 x 10 15 f ~rzr~7~~77 ) (1 + z)^ 2 M Q (19) 



The corresponding scaled total mass profile, M to t(x), is 
directly related to the variation of the gas density slope 
a n (x): 

M tot (x) 



M to t(x) 



M 



VT200 



0.472 x a n (x) 



(20) 



Fig. [j] shows the total mass profiles, M tot (x), derived from 
our best fit /3-model A (dotted lines) and B (full lines), 
for all the clusters of the spectroscopic subsample. As for 
a n (x), the differences between the two models are smaller 
than the dispersion in the cluster population. 

The scatter in the total mass profiles is the same as 
for a n (x) and is minimal at about x s = 0.3, the specific 
radius of the /3-core radius relation. At that scaled radius 



3 assuming that the cluster is in a dynamical relaxed state 
and approximately spherically symmetric 
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Fig. 10. Comparison between mass estimates for the spec- 
troscopic subsample. Bold line: Numerical simulations 
of EMN96. Full lines: The middle curve is the mean 
mass profile for our sample derived from the isothermal 
/3-model (x c = 0.13 and (3 — 0.67); The top and bot- 
tom curves correspond to this mean plus or minus the bi- 
weight dispersion. Dotted line: The mean mass profile for 
the same /3-model but with a polytropic temperature pro- 
file of index 7 = 1.2, approximating the results of Marke- 
vitch et al. (1998). 



a a (x s ) = 1.67 ±0.18 and M to t(x s ) = 0.236 ± 0.025 
physical coordinates this corresponds to: 

1/2 

r,. = — I • : I Mpc 



Mtot(r s ) 



(l + z)3/ 2 VWkeV 
(6.90 ±0.74) x 10 14 



In 



(21) 



3/2 



Mr. 



(22) 



(1 + z) 3 / 2 V 10keV , 

Note that the numerical values in these equations do not 
depend of the calibration of ryT200 & n d Mvt200 ■ The ra- 
dius of minimal dispersion corresponds to a density con- 
trast of S = 1750 (for the mass computed from the isother- 
mal /3-model ). 

The mean mass profile, corresponding to x c = 0.13 
and (3 — 0.67, and the two curves corresponding to this 
mean plus or minus the bi-weight dispersion are plotted 
on Fig [H]. For comparison we also plotted the mass profile 
derived from the numerical simulations of Evrard, Met- 
zler & Navarro (1996; hereafter EMN96; Table 5). The 
isothermal /3-model yields lower values but the discrep- 
ancy decreases with radius and is less than 10% around 
x — 1. Both models would give equal values at x = 1 for 
/3 ~ a n (l)/3 = 0.71 (from Eq. ^0|), whereas the mean /3 
value is only 6% lower. 

The isothermal assumption may give poor estimates of 
the total mass if strong temperature gradients are present. 
Recent observations with ASCA suggest that the tempera- 
ture does decrease with radius, although the uncertainties 
on the temperature profiles are still large (Markevitch et 



al. 1998). To illustrate the impact of such a temperature 
gradient we considered the composite projected temper- 
ature profile derived by Markevitch et al. (1998, Figure 
7) from a sample of bright clusters. This profile is steeper 
than most of the profiles derived from available simula- 
tions (Markevitch et al. 1998) and thus is likely to give an 
upper limit on the influence of a temperature gradient, at 
least up to x — 0.5 the limit of the observations. We found 
that the shape of the projected temperature profile, T{x), 
can be well approximated by a polytropic model with an 
index of 7 = 1.25, a normalisation factor A — 1.2 and our 
mean /3-model density profile (/3 = 0.67, x c = 0.13): 

Jw.^ir'-Afi+w*)'!-*' 7 -" (23 , 



= A 



Tx V n go, 
Assuming that the observed temperatures correspond to 
the emission weighted mean temperatures along the line- 
of-sight, it is easy to show that the radial temperature 
profile, T(x), also obeys the same polytropic law with a 
higher normalisation factor, about 10% higher for the pa- 
rameters we derived. 

For a polytropic gas, the temperature logarithmic slope 
ax, entering the hydrostatic equation (Eq. 17), is simply 
related to the corresponding density slope: 
a T (x) = (7-I) a n (x) (24) 
The corresponding mass profile, plotted on Fig. [H], is thus: 
T(x) 



M tot (x) = 0.472 7: 



x a n (x) 



(25) 



The mass is larger than the isothermal estimate be- 
low a; ~ 0.35 and higher above. At x = 0.3 the correction 
is only 4%, less than the typical scatter. At large radii 
the differences are drastic: at x = 1 the isothermal esti- 
mate is 1.7 times higher than the estimate using the poly- 
tropic model. Whereas in the isothermal case we deduce a 
mass profile shallower than found in numerical simulation 
(EMN96), the mass profile is now slightly more peaked 
than the theoretical one. We must however emphasize that 
the temperature profile used is based on an extrapolation 
of the data beyond x = 0.5. 



8. Discussion and Conclusions 

We have studied the shape of the surface brightness pro- 
files of a sample of clusters covering a large variety of 
morphological types. The surface brightness distribution 
in the ROSAT band hardly depends on temperature and 
can be directly translated into the gas density distribution. 
Although the gas density profiles of the clusters are not 
perfectly similar, our quantitative analysis gives evidence 
for structural regularity: 

— We found converging evidence that there is a scaling 
radius in clusters which varies as T 1 / 2 . The surface 
brightness profiles appear more similar once the phys- 
ical radii are scaled. The /3-core radius relation is tight- 
ened and the scatter of the density profile slopes is re- 
duced. However the small range in temperature and 
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rcdshift of our sample precludes any definitive conclu- 
sion, in particular on the exact variation of the scale 
radius with T. 

— There is a parabolic correlation between the two shape 
parameters of the /3-model : the core radius and the 
slope. Cluster density profiles essentially constitute a 
one shape parameter family. 

— We observe a very large scatter of the surface bright- 
ness profiles in the central part of clusters. A central 
excess compared to the /3-model can be present, typi- 
cally within x C ut = 0.08. In this region the logarithmic 
slope can vary by more than a factor of 2 and the emis- 
sion measures at the center are spread over two orders 
of magnitude. 

— Beyond a scaled radius of x = 0.2, the density pro- 
files are remarkably similar in shape and resemble a 
/3-model profile with a core radius of x c = 0.13 and 
= 0.67. The bi-weight dispersion of the logarithmic 
slope a n (x) at a given radius is less than 20%. At the 
virial radius the logarithmic slope is close to 3/3. Its 
distribution reflects the distribution of this parame- 
ter, which has a bi-weight dispersion of only 18% but 
presents a tail of rare high values. The scatter is min- 
imal at x = 0.3 where a n (x) — 1.67 ± 0.18 with a 
gaussian distribution. This minimum is related to the 
correlation between x c and 0. These results at x = 0.3 
are insensitive to the exact treatment of a possible cen- 
tral excess when fitting the profiles. 

As the gas evolves in the potential of the dark matter 
it is unlikely that such regularity in the gas density profile 
shape can be observed without a comparable regularity 
in the dark matter profiles. Similarity of the dark mat- 
ter profiles is naturally expected in hierarchical clustering 
scenario (Teyssier, Chieze & Alimi 1997; Navarro, Frenk 
& White, 1997) and our results therefore give support to 
it. However, the gas density is not the sole tracer of the 
gravitational potential. If a universal profile for the dark 
matter is responsible for the similarity of the X-ray pro- 
files, one also expects that the temperature profiles have 
similar shapes. Both numerical simulations (e.g. Navarro, 
Frenk & White 1995; EMN96) and recent X-ray observa- 
tions (Markevitch et al. 1998) indicate that this is proba- 
bly the case. 

Our results are consistent with a scaling radius which 
goes as T 1 / 2 as expected if clusters form at constant over- 
density. However our sample is too small to exclude other 
scaling laws as, for instance, the slight variation of over- 
density with mass obtained by Navarro, Frenk & White 
(1997). 

Similarity breaking in the X-ray properties of clusters 
is expected if non adiabatic processes play an important 
role in the evolution of the gas. These include radiative 
cooling and possible heating by galactic winds. 

Radiative cooling is known to be more important in 
the cluster core where the cooling time, which scales as 



the gas density, is shorter. The presence of Cooling Flows 
of various strengths most probably explains the variation 
we observe in the density profiles below x = 0.1. The 
large scatter could be either due to intrinsic variation in 
the initial gas density in the core, and thus in the cool- 
ing time, or to variation in the Cooling Flow "age" . The 
first explanation is unlikely in view of the observed quasi- 
similarity of the scaled emission measure profiles at large 
radii. Let consider perfectly similar clusters obeying the 
assumptions used in the scaling of the profiles. The mean 
density of clusters in a given redshift range is a constant 
(formation at fixed overdensity) and since the gas mass 
fraction does not vary from cluster to cluster, the mean 
gas density is also a constant. As the density profiles are 
similar in shape, the central density is also constant and 
so is the cooling time. No major break of similarity in the 
core compared to external regions is thus expected, if clus- 
ters have the same "age" . The very large scatter observed 
in the core properties thus favors a scenario where Cool- 
ing Flows are recurrent phenomena that are periodically 
erased by strong mergers, a natural feature of hierarchical 
clustering (Fabian et al. 1994). In that case the observed 
variety of core properties in the cluster population would 
naturally reflect the statistics of the formation process via 
merger events. 

Non gravitational heating, such as extra energy inputs 
by galactic winds, is thought to play an important role in 
the ICM physics. In " pre- heating" scenarios, early winds 
add roughly a constant entropy to the system and thus 
affect cooler clusters more strongly, hence breaking the 
similarity law. The main expected effect of these winds is 
to inflate the gas distribution. Preheating scenarios can 
thus explain the observed steepening of the Lx~T and 
size-temperature relations relative to the similarity scaling 
(Arnaud & Evrard 1998, Mohr & Evrard 1997, Cavaliere, 
Mcnci & Tozzi 1998). Further support to this scenario was 
recently given by Ponman, Cannon & Navarro (1998) who 
showed that the entropy of cool systems is indeed higher 
than what can be achieved through gravity alone. They 
consistently put into evidence similarity breaking in the 
surface brightness profiles of clusters: cooler clusters have 
systematically shallower profiles than hotter systems. The 
regularity we found in the gas density profiles beyond the 
core does not contradict such results. Our spectroscopic 
subsample comprises relatively hot clusters. Our quanti- 
tative study simply shows that non gravitational heating 
is negligible above kT ~ 3.5 keV. This is in agreement 
with the quantitative deviations from similarity laws ob- 
served in clusters. The gas entropy at x = 0.1 plotted 
by Ponman, Cannon & Navarro (1998) as a function of 
cluster temperature levels off from the similarity law only 
below ~ 3 keV. Similarly the variation of the structure pa- 
rameter Q{T) = (n gas ) I (n gas ) 2 , which characterises the 
concentration of the gas distribution and plays a role in 
the steepening of the L^-T relation, has been studied by 
Arnaud & Evrard (1998, Figure 2). An increase with tern- 
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perature has been noted but the effect is again stronger 
below 4 keV. The existence of such a "threshold" around 
3 — 4 keV clearly favors scenarios with constant entropy 
inputs. In turn it shows that the outer regions of hot clus- 
ters are good tracers of purely gravitational processes. 

The mean gas density profile we derived is shallower 
than found in CDM simulations without winds. For in- 
stance Navarro, Frenk & White (1995, Table 2) found a 
mean core radius x c = 0.10 and a mean /3 of 0.82, while 
Eke, Navarro & Frenk (1998, Table 3) derived for a low 
VI Universe average values of x c = 0.05 and /3 = 0.73. If 
pre-heating is negligible in our spectroscopic subsamplc 
other effects can explain the difference. First the transfer 
of entropy from the dark matter to the gas, which has been 
invoked to explain the segregation between these two com- 
ponents (Eke, Navarro & Frenk 1998; see also Teyssier, 
Chieze & Alimi 1997), may be larger than predicted in 
the simulation. We can also speculate that the dark mat- 
ter distribution itself is actually less concentrated, as, for 
instance, expected if some hot dark matter is present (see 
also below). 

The gas structural variations throughout the cluster 
population are small beyond x = 0.2 but depend on the 
radius considered. The logarithmic slope a n (x) of the den- 
sity profile presents a minimum spread at x ~ 0.3. This 
scatter then increases with radius and some clusters present 
large asymptotic slopes compared to the mean asymptotic 
slope. The /3-model , used in our study, is reasonably accu- 
rate to estimate both the gas density distribution and total 
mass distribution, as shown by the hydrodynamic simula- 
tions of Schindler (1996) or EMN96. However, by fitting 
such a model to simulated galaxy clusters, these last au- 
thors found that it introduces an extra-scatter in the total 
mass estimate and that this scatter increases with radius. 
We probably see the same effect here. In the hydrostatic 
isothermal /3-model , the logarithmic slope of the density 
profile is directly linked to the total mass. For a given dark 
matter profile, scatter in the gas density slope estimate can 
be introduced by i) errors in the measure due to the back- 
ground determination and statistical errors; ii) departure 
from a spherically symmetric /3-model for the gas distri- 
bution due to imperfection of the model, bimodality, an- 
cient mergers, and high eccentricities; iii) departure from 
the hydrostatic equilibrium. The first and third effects are 
expected to be larger in the outer regions. On the con- 
trary, the signal to noise is very high in the region around 
x = 0.3 and the gas density slope, which is the derivative 
of the profile at that radius, is tightly constrained. This 
probably explains why the scatter in the slope estimate is 
minimum there and why the derived slope value is insen- 
sitive to the exact model used to fit the data (model A 
or B). In conclusion the observed variation of the scatter 
with radius is more likely an artifact due to the modeling 



of the data, rather than the direct consequence of similar 
variations in the dark matter profiles. 

The slope at x ~ 0.3 is almost the same for all clus- 
ters. Since, in the /3-model , it depends on both (3 and 
core radius, the values of these two parameters are cor- 
related throughout the cluster population, as we did ob- 
serve. This very correlation between /3, that controls the 
outer slope, and x c , that controls the core size, may in- 
dicate that the best model to fit the data is not the two 
parameter /3-model but some other model, involving only 
one shape parameter (scaling with the virial radius). 

As a result of the small scatter on a n (x), the spread 
in the scaled mass profiles derived from the hydrostatic 
isothermal /3-model is small. The uncertainty on the dark 
mass profile, deduced from the hydrostatic equation, is 
dominated by the systematic uncertainty on the temper- 
ature gradients, rather than by the intrinsic scatter in the 
gas density profiles. 

The radius at which the scatter of the gas logarithmic 
slope is lowest is most likely in the region in which the 
/3-model gives the most accurate result on the total mass 
of clusters. It is also in the region in which temperature 
gradients are not expected to play a large role, as indicated 
both by numerical simulations (e.g. EMN96) and by the 
observed temperature profiles (see previous section). We 
thus propose the M-T (r-T) relation, we derived at x = 
0.3 (Eq. p^) from the hydrostatic isothermal /3-model , as 
a reference test point for numerical simulations. 

In this relation the mass scales classically as T 3 / 2 (1 + 
z) -3 / 2 . The normalisation at x = 0.3 is however ~ 30% 
lower than found in the numerical simulations of EMN96, 
whereas these authors did not observe any significant bias 
at that radius in the isothermal /3-model estimate. One 
can speculate again that the dark matter profile is in 
fact shallower than in CDM simulations, as expected if 
some HDM is present. However the discrepancy could also 
indicate that the temperature spatial variations or bulk 
motions and residual velocity dispersion of the gas are 
higher than expected. In that case the hydrostatic isother- 
mal equation underestimates the real mass. Large velocity 
fields are naturally associated with recent formation or ac- 
cretion. Therefore their properties should strongly depend 
on the specific history of each individual cluster and the 
corresponding velocity profiles are expected to present a 
high relative scatter. If large kinetic energy is common in 
the cluster population, it would thus be surprising that the 
gas density profiles remain similar. In conclusion we think 
that the regularity of the gas density profiles supports the 
validity of the hydrostatic equation. On the other hand 
quasi-similar temperature profiles corresponding to large 
temperature gradients could well be present. It thus is es- 
sential to accurately measure the temperature profiles of 
large sample of clusters to establish firmly the normalisa- 
tion of the M-T relation and the shape of the dark matter 
profile. This will be soon possible with XMM and AXAF. 
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